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EXPERIMENTAL AND ANALYTICAL INVESTIGATION OF INTERFACIAL 
HEAT AND MASS TRANSFER IN A PRESSURIZED TANK 
CONTAINING LIQUID HYDROGEN 
by William A. Olsen 
Lewis Research Center 

SUMMARY 

Temperature measurements in the vicinity of the placid interface of two-phase hy- 
drogen contained in a pressurized tank have yielded information on interfacial heat and 
mass transfer. These measurements indicated that interfacial mass transfer can have 
an appreciable effect on bulk gas and liquid conditions, while heat transfer from the gas 
can probably be neglected in most cases. An analysis of interfacial heat and mass trans- 
fer, which closely predicts the data, is included herein. 


INTRODUCTION 

Cryogenic propellant tanks are pressurized to cause outflow, prevent pump cavita- 
tion, and/or stiffen the structure of the tank. When a tank is suddenly pressurized 
(artificial pressurization) to a new level from an initial equilibrium state, the two-phase 
fluid system attempts to go to a new equilibrium situation at the higher pressure by the 
transfer of heat across the tank walls from the environment, and heat and mass transfer 
across the interface. Heat transfer from the environment heats liquid at the tank walls, 
which then flows up the walls and forms a growing layer of relatively warm liquid below 
the interface. Heat and mass transfer across the interface tends to enlarge this warm 
layer of liquid. Similar energy and mass transfer occur while the liquid is expelled 
from the tank or the tank pressure builds up slowly as a result of the heat leak (self- 
pressurization). 

There have been several analytical and/or experimental studies that have examined 
this complex nonequilibrium process (refs. 1 to 10). However, none of these studies 
experimentally investigated the interfacial region and interfacial heat and mass transfer 
in detail. 



In references 1 and 2, analytical and experimental studies of pressurization were 
performed. The former used an overall heat and mass balance that yielded only rough 
estimates of interfacial heat and mass transfer. References 3 and 4 were concerned with 
stratification in high-heat-flux tanks but did not consider interfacial transfers. In refer- 
ence 5 the effect of interfacial transfers on tank pressure was experimentally observed 
when the bulk liquid or gas was artificially mixed. Similar analyses of interfacial heat 
and mass transfer were performed in references 6 to 8. In reference 9, pressurization 
and stratification were studied analytically considering interfacial transfers. Probably 
the most applicable experimental study of interfacial transfer is given in reference 10, 
which obtained detailed experimental temperature profile data T = T(Z, t) for a very low- 
heat-flux self-pressurized tank. (All symbols are defined in the appendix. ) Only a very 
limited evaluation of interfacial heat and mass transfer was performed. 

As part of a general study at the Lewis Research Center of internal fluid conditions 
in pressurized propellant tanks containing cryogens, the interfacial region was experi- 
mentally studied more extensively than in previous references to determine the inter- 
facial heat and mass transfer for artificially and self-pressurized tanks where the heat 
flux was not small. In addition, the results of this experimental study were compared 
with corresponding results of the now standard analytical method to predict heat and mass 
transfer in suddenly pressurized tanks. 


ANALYSIS OF INTERFACIAL HEAT AND MASS TRANSFER 
Physical Description of Problem 

When a vented tank containing a two-phase one- component cryogenic fluid in equilib- 
rium is rapidly pressurized to a new constant pressure by adding gas, the system tends 
toward a new equilibrium condition at that new pressure (fig. 1). Gas near the interface 
is rapidly compressed by the well-diffused (little mixing) incoming gas. Gas near the 
interface, which was heated by compression, begins to cool by heat transfer to the liquid. 
The interface immediately assumes the new equilibrium state (saturation). (Kinetic 
theory requires some difference between the interface temperature and saturation; how- 
ever, this difference is generally small for a pure fluid, refs. 6 to 8. ) Concurrently, 
the liquid begins to be heated by heat and mass transfer from the warm gas and heat 
transfer from the environment. 

All the heat transferred to the liquid will be stored nonuniformly in the liquid, re- 
sulting in stratification (i. e. , temperature gradients, T = T(Z, t) ). The environmental 
heat flux through the side walls of the tank to the liquid causes relatively warmer liquid 
to flow up the tank and form a warm layer near the interface (fig. 1), which is further 
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heated by interfacial heat and mass transfer. Heat transferred across the tank bottom 
from the environment will essentially heat up the bulk liquid below this warm layer of 
liquid (ref. 11). 


Analysis of Artificial Pressurization 

Except for the specification of some of the boundary conditions, the analysis of this 
problem has been performed many times (refs. 6 to 8). The analysis used in reference 6, 
which is most easily used here, is summarized subsequently only for the reader’s con- 
venience. Basically these analyses assume that the flows and heat transfer occurring in 
the region of the interface of the complex system previously described can be physically 
approximated by two semi-infinite volumes of gas and liquid joined by an interface (fig. 2). 
Using a moving coordinate system, attached to the interface, gives the differential equa- 
tions and boundary conditions describing conditions in the gas and liquid near the inter- 
face as 
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The boundary conditions are (fig. 2): 




t < 0, 

Z = Z: 

T f = T g = T o 

(3), (4) 

t> 0: 

'"Z = +°o ; 

T * =T b 

(5) 

< 

Z = -°°: 

T g = T r 

(6) 


jz = 0: 

T g = T , = T i 

(7), (8) 


where T q , T^, T f , and are not functions of time or position. 

The velocities and are obtained by an interfacial mass balance and an 
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interfacial energy balance. The interfacial mass balance, which satisfies the continuity 
equation is 


W 



The interfacial energy balance is 
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( 9 ) 


( 10 ) 


The partial differential equations (eqs. (1) and (2)) are transformed to ordinary dif- 
ferential equations by the following substitutions (refs. 6 to 8) 



: Pj< T i - T t> 


(id 



and the similarity substitution 



( 12 ) 


(13) 


where j denotes gas or liquid. These substitutions (eqs. (11) to (13)) convert the par- 
tial differential equations (1) and (2) to ordinary differential equations: 
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The resulting transformed boundary conditions become: 
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(16) 

II 

8 

<£> 

aq 

1! 

(17) 

'g =0: e ( = 0 g =o 

(18), (19) 


The interfacial energy and mass balances are transformed by equations (11) to (13) to 


p 4 (« £ ) 1/2 v £ = P g <V 1/2 v g 


( 20 ) 
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Solutions of the differential equations (eqs. (14) and (15)) are given subsequently. The 
solutions describe the temperature profiles in the gas and liquid in the vicinity of the 
interface: 


e = e 

g r 


e * = *b 


erf(V g ) - erf(V g - ri g ) 
erf(V g ) - 1 


erf(V £ ) - erf(V^ - 7] £ )" 


erf^) + 1 

Interfacial mass transfer is obtained from equations (9) and (12): 



( 22 ) 


(23) 


(24) 


Interfacial heat transfer to the liquid is obtained from equations (11), (13), and (23) ac- 
cording to 


5 



to achieve 



Interfacial heat transfer from the gas is obtained from equations (11), (13), 
cording to 


to obtain 
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The nondimensional velocities Vg and required in equations (24), (26), 
obtained from the simultaneous solution of 


-7T 


1/2 h 


(v t [e X p(v f 2 )]}[l + erf(V,)] {v g [exp(v|)]}[erf(V g) - l] 


and 


P f <a ( ) 1/2 V, = P g ^/ /2v g 


( 26 ) 

and (22) ac- 
(27) 


(28) 

and (28) are 
(29) 


( 20 ) 


Equation (29) was obtained by differentiating equations (22) and (23) and substituting into 
equation (21). 



Values of the constant thermal boundary conditions (eqs. (5) to (8)) will now be spe- 
cified based on the physics of the problem. The interface temperature should be 
very nearly the saturation temperature =* T ga j., which will be constant after the tank 
has been rapidly pressurized to a constant pressure. The interface is assumed to offer 
no resistance to mass transfer (refs. 8 and 12). 

The initial gas is assumed to be rapidly compressed by the incoming gas with no 
mixing, and thereafter the heat transfer is small so that the changes in T r are small. 
A reasonable choice for T , therefore, appears to be the adiabatic compression tem- 
perature of the gas, which will be constant. Thus T r is assumed to be 


ad 



(30) 


where y ^ 1. 8 for gaseous hydrogen for 37° R < T r < 50° R. The bulk liquid boundary 

condition T^ is obtained by assuming that the change in the liquid temperature through 

the heat leak is small enough so that T^ can be taken as the initial temperature T . 

Thus, T k ^ T . 

’ b o 

In summary, the values for the boundary conditions are considered to be 


^i ^sat ^sat(^ 


(31) 


y-1 


T r = T ad = AdW = t o(^T r 

Vj V 


T b = T o 


(32) 

(33) 


and the properties can be evaluated at T^ = T ga ^(P) since the analysis only holds near 
the interface. 
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APPARATUS AND PROCEDURE 


General Description of Experiment 

The experiment can be described as a tank containing a two-phase fluid initially at a 
uniform temperature and low pressure that is suddenly raised to a higher pressure. The 
pressure is either increased quickly by the addition of gaseous hydrogen or helium 
(artificial pressurization) or slowly by heat leak (self-pressurization). Fluid tempera- 
tures, tank pressure, and liquid level data are recorded as the fluid tends toward a new 
equilibrium state. Tank pressure, pressurant gas, and method of pressurization are 
varied as parameters (see table I). 


Description of Test Facility 

The test facility consisted of a large vacuum- jacketed tank (30 in. in diam. by 88 in. 
long, fig. 3(a)) that could be filled, emptied, and pressurized automatically. A 3 -foot- 
long carbon- resistor temperature rake was vertically fastened to the lid so that it was at 
the center of the tank. A viewing window, allowed the experimenter to observe the fluid 
conditions in the tank. The tank could be pressurized up to 40 pounds per square inch 
gage. 


Description of Temperature Rake 

The 3 -foot-long vertical temperature rake used in this study consisted of forty-two 
0. 07 -inch- diameter 1/10 watt calibrated carbon resistors, which were fastened to a 
Bakelite superstructure at various levels (fig. 3). Thirteen resistors, which supplied the 
significant interfacial data, were spaced 1/4 inch apart at the middle of the rake (fig. 3). 
An effort was made to keep the carbon resistors isolated from significant thermal and 
hydrodynamic error sources by means of low conductivity supports and lead wires (man- 
ganin alloy) and a tree-like rake structure that did not disturb the flow significantly. 


Procedure 

Two types of tests were performed. One was artificial pressurization, where gas 
at about room temperature is forced through a diffuser into a tank to bring the tank to 
some higher pressure, which is then held constant (tests 1 to 3). The diffuser was de- 
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signed so that it caused very little gas mixing. The other type of test was self- 
pressurization, where the vent valve is closed and the tank pressure would build up by 
heat leak (test 4). 

Artifi cial pressurization (tests 1 to 3) . - The liquid level was initially visually set at 
about the middle of the 1/4 inch spaced resistors on the temperature rake. The initially 
vented tank was then quickly pressurized with the desired gas to a preset pressure, which 
was thereafter held constant. About 1 minute was required to attain the desired pressure. 
At specific time intervals after the start of pressurization, the temperature profile and 
tank pressure were recorded by automatic self-balancing digital recording equipment. 
After completion of a test, the tank pressure and/or pressurant gas were then changed, 
and the preceding procedure was repeated. Throughout each test the liquid condition was 
visually monitored. The heat flux was measured for a given liquid depth and vacuum 
pressure by noting the boiloff before and after each run. 

Self-pressurization . - The only difference between the self-pressurization and the 
artifical pressurization tests was that the vent valve was closed to build up pressure 
slowly by heat leak in the former whereas pressurant gas was added in the latter to at- 
tain a constant pressure. 


Instrumentation Accuracy 

The position of the carbon resistors with respect to one another was known to within 
0. 05 inch including thermal dimension changes. The entire temperature rake was im- 
mersed in boiling liquid hydrogen (at atmospheric pressure) before each test run to 
check its thermal accuracy at a known condition. Temperature errors due to thermal 
conduction through the leads and the rake, thermal radiation from the tank walls and lid, 
and electrical heating in the carbon resistors were estimated analytically. Based on 
these estimates and the temperature check in boiling hydrogen, it is expected that a ther- 
mal error of no greater than ±0. 1° R will result. The thermal time lag (time constant) 
of the carbon resistor thermometers in liquid and gaseous hydrogen (maximum tempera- 
ture difference < 200° R) was analytically estimated to be less than 0. 2 and 2 seconds, 
respectively. This lag results in a negligible thermal error for this slow transient ex- 
periment. 

The automatic system that controlled the tank pressure maintained a fairly constant 
pressure that varied no more than ±5 percent around the nominal pressures listed in 
table I. Tank pressure was sensed within ±0. 1 percent of full scale, which would result 
in an error in the saturation temperature no greater than ±0. 05° R. 

The overall heat flux (based on the liquid- wetted area) for a given liquid depth and 
vacuum jacket pressure was determined by measuring the rate at which liquid was 
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boiled off when the tank was vented to the atmosphere (P' = 0 psig). From this informa- 
tion the heat flux could be approximately calculated by using equation (34): 



(34) 


Because of greater possible changes in the liquid temperatures at the higher tank pres- 
sure levels, it is recognized that the experimentally determined heat-leak values will be 
in excess of the actual quantities (i. e. , (Q/A) w is proportional to - T^ ). Within the 
limits of the tests conducted, however, this error should not exceed 2 percent. 


DATA PROCESSING 


The interface location is determined from temperature profile data (T = T(h, t) ) by 
assuming that the interface is saturated. This assumption was also made in the analyses 
and in other references (e. g. , refs. 1, 6 to 8, and 10). From the level location and tem- 
perature profile data, interfacial data can be obtained. For example, the interfacial heat 
transfer from the gas is obtained from equation (25) where (<lg/A) int = - (q^/A^ and 
STg/aZ = -(3Tg/ah) to give " " b 



(35) 


The interfacial heat transfer to the liquid is similarly obtained from equation (27) to give 



The mass flux transferred across the interface for a given time can be obtained ex- 
perimentally from the one-dimensional energy balance described by equation (10) where 



and 
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It is assumed that 


which gives 



(37) 


Two further assumptions are inherent in equations (35) to (37). Radial temperature 
gradients in the gas and liquid in the vicinity of the interface (when the region near the 
tank walls is excluded from consideration entirely) were assumed to be negligible. Ra- 
dial temperature profiles near the interface have been measured in pressurized outflow 
studies (unpublished data obtained by R. DeWitt of Lewis). These data indicate that the 
assumption is reasonable. The second assumption is that the significant temperature 
gradients occur over a sufficient distance so that a finite number of thermometers 
(carbon resistors spaced 1/4 in. apart) would be able to detect and to measure these 
gradients. 


RESULTS AND DISCUSSION 

Vertical temperature profiles in the liquid and gas and tank pressures were taken at 
various time intervals after the start of pressurization. The method of pressurization, 
pressurant gas, and pressure level were varied (see table I). 


Artificial Pressurization with Hydrogen Gas Tests - 1 and 2 

Raw data . - In the interest of brevity, only one typical group of rake-length tempera- 
ture profiles (T = T<h, t, P'>) will be reported. Figure 4 contains the profiles for test 1 
at a tank pressure of 40 pounds per square inch gage. The portion of this profile where 
the temperature sensors are spaced at 1/4-inch intervals (i. e. , the region about 1. 5 in. 



on either side of the interface) is enlarged in figure 5(c). Figures 5(a) and (b) contain 
this type of data for tank pressures of 8 and 20 pounds per square inch gage. Figure 6 
contains a similar enlargement for test 2 data. The solid lines represent the author’s 
judgement of the best fairing. The saturation-dictated interface location is also plotted 
in each of these figures. 

Interfacial heat and mass transfer. - Interfacial mass transfer, as obtained from the 
analysis* (eqs. (20), (24), and (29)), is plotted in figure 7. The experimental interfacial 

mass fluxes that occurred during tests 1 and 2 are determined from the temperature pro- 
files (figs. 5 and 6) at the saturation-dictated interface according to equation (37). These 
data are also plotted for comparison in figure 7. Agreement between the analysis and 
the experimental data is generally good. Figure 7 indicates that the mass transfer is 
condensation which increases with tank pressure and decreases with time. 

The analytically determined interfacial heat transfer to the liquid (eqs. (20), (21), 
and (26)) is plotted in figure 8. Also shown are experimental values of heat transfer to 
the liquid as determined from equation (36). These results also indicate good agreement 
between the analysis and the experiment. Again, interfacial heat transfer to the liquid 
increases with pressure and decreases with time. 

Figure 9 contains the analytical (eqs. (20), (21), and (28)), and experimentally 
(eq. (35)) determined values for interfacial heat transfer from the gas when it is assumed 
that there is no discontinuity in the slope at the interface i. e. , (9Tg/9Z)^ = (9Tg/ 9Z)^ .J 
Agreement between the analysis and the experiment is poor especially for increasing tank 
pressure. Fortunately, heat transfer from the gas is small compared with heat transfer 
to the liquid so that mass transfer is not greatly affected by this error. This can be seen 
from equation (37) where kg » k^. The reasons for this disagreement are discussed 
in the section titled Interfacial boundary conditions. 


Artificial Pressurization with Helium - Test 3 

In a two- component two-phase system, interfacial heat and mass transfer caused by 
temperature differences and component mass transfer (diffusion) caused primarily by 
concentration differences are coupled. Consequently, it should be expected that experi- 
mental measurements would be complex. Because of the scarcity of diffusion- coefficient 
data and the difficulty in sampling the concentration of components in a transient system, 
only significant raw data from the helium pressurant test (test 3) will be reported. The 

*A copy of the analytical computer program, written by Margaret M. Smith, can be 
obtained from Instrument and Computing Division, Lewis Research Center. 
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temperature profiles near the interface are plotted in figure 10. 

In general, the effective saturation temperature of this two-component system 
(hydrogen- helium pressurant) T gat * should not be the same as the saturation temperature 
for the one- component system (hydrogen, tests 1 and 2) T g ^. As a rough indication of 
the interface location in figure 10, it was assumed that T ga ^. T ga j.*; however, no more 
can be done with this data because of the lack of concentration data. 


Self- Pressurization - Test 4 

With the vent valve closed, the tank pressure rose from P* = 0 pounds per square 
inch gage at a uniform rate of 0. 4 pound per square inch gage per minute as a result of 
the heat leak. Figure 11 is a plot of the temperature profiles in the vicinity of the inter- 
face for this test. A comparison of heat and mass flux data between artificial and self- 
pressurization tests reveal different trends; in this test, both quantities increased with 
time. Some slight boiling was noted at the tank walls near the interface during most of 
the test. Excluding the mass transferred by this boiling at the walls, heat and mass 
transfer was calculated by equations (35) to (37) and plotted in figure 12. These fluxes 
are much smaller than those for artificial pressurization. The analysis is not applicable 
for this case where the saturation temperature varies with time because of the time- 
varying pressure, therefore no effort was made to correlate these data with the analysis. 


Discussion of Boundary Conditions 


The boundary conditions were specified on the basis of physical arguments. The ob- 
ject of this section is to discuss and justify, by means of analysis and experiment, the 
values of T r and T^ and the interfacial boundary conditions. 

Bulk fluid boundary conditions . - The bulk liquid temperature, which was assumed 
to be constant (i. e. , T^ = T ) in the analysis, will actually increase slowly with time 
because of the heat leak from the environment to the liquid. If it is assumed that all the 
heat leak from the environment is stored uniformly in the bulk liquid, T^ can be approxi- 
mated by 


T b = “ 






(38) 
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This function (eq. (38)) is approximated by a series of finite steps and is used in the 
analysis to give an estimate of the effect of this time-varying boundary condition on inter- 
facial mass transfer. The result of this calculation is plotted in figure 13 for pressures 
of 20 and 100 pounds per square inch gage and various values of 

k|k= (Q/A) w /3600p £ Cp d] and H^/D. The two single-line curves, 5 and 1 are for tanks 

with zero environmental heat flux (K = 0). The band curves 2 to 4 and 6 to 8 describe the 
range of mass fluxes for 2 > H^/D > 0. 5 that occur for various heat fluxes. The upper 
part of the band corresponds to an H^/D of 2 and the lower part to an H^/D of 0. 5. 

Figure 13 indicates that the effect of environmental heat leak on interfacial mass 
transfer is to decrease the mass flux to a lower value than would be predicted by a con- 
stant value of T^(T^ = T Q ). Compare curve 5 to 6 and 7, or compare 1 to 2 and 3 in fig- 
ure 13. This difference increases with time. The effect on interfacial mass transfer by 

-4 o 

changes in T K caused by heat leak is small for t < 10 minutes and K < 7x10 R per 
D 5 

second, or t < 30 minutes and K < 7x10 (fig. 13). The former situation would cor- 
respond to a heat flux of 85 Btu per square foot per hour for a 3 -foot-diameter tank or 
about 280 Btu per square foot per hour for a 10-foot-diameter tank. The M centerline” 
of the band of curve 2 corresponds to the experimental situation of tests 1 and 2 where the 
heat flux was about 8 Btu per square foot per hour. It appears that interfacial mass 
transfer can be predicted approximately by using the constant liquid boundary condition 
(T^ = T ) for most practical situations (where K = (Q/A) 
and the time is small). 

Some justification is required for the choice of the bulk gas thermal boundary condi- 
tion T r chosen in the analysis. Both the analysis and experiment were based on situa- 
tions where the gas volume is not appreciably mixed by the rapid introduction of pres- 
surizing gas. In this situation, one might expect the pressurizing gas to compress the 
original gas adiabatically during a rapid pressurization. The main justification for the 
choice of the boundary condition T f = T ac j is that the interfacial heat transfer to the 
liquid and mass transfer predicted by the analysis closely agrees with the experiment. 

On the negative side there is poor agreement in interfacial heat transfer from the gas. 
This disagreement is discussed further in the next section. 

Interfacial boundary condi tions. - Based on the usual thermodynamic assumption of 
local equilibrium, the following are the accepted theoretical interfacial boundary condi- 
tions (ref. 13): 


w /3600 p 


f c D <7X 10 - 


°B/sec 
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+ p W L = 0 
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References 6 to 8 made use of these two equations in their analysis, and they have been 
used in this report in order to evaluate the data. The first equation (eq. (7)) is considered 
approximate and should be experimentally verified. Considering the temperature profile 
data it seems reasonable that (0Tg/9Z^ ~ (dT dzj^, and since is about an order of 
magnitude greater than kg, it is reasonable to state that heat transfer from the gas 
(first term in eq. (10)) has little effect on mass transfer (third term) or on heat transfer 
to the liquid (second term). Therefore it should be expected that the data analysis, which 
assumed ^9Tg/9zj ^ = (dT^/dZ^, should give good results for interfacial mass transfer 
and heat transfer to the liquid, but poor results for interfacial heat transfer from the gas. 
This was indeed shown to be the case in figures 7 to 9. The analysis indicates that 
(aTg/azji s (aiyaz). is not a good assumption. Figure 14 is a plot of the analytically 
determined ratio of the liquid and gas temperature gradients at the saturation-dictated 
interface (aTg/az) i/(pT^/9Z^. Such small discontinuities could not have been ob- 
served in this experiment if they had existed, because the sensors are 0. 07 inch in diam- 
eter and 1/4 inch apart. Examples of such local slope discontinuities are plotted in fig- 
ure 5. 


Effect of Interfacial Heat and Mass Transfer on Bulk Fluid Conditions 


Although the flow of condensate appears to be small (fig. 7), its effect on conditions 
in the bulk gas and liquid may not be small in the artificially pressurized case. For ex- 
ample, consider a gas space (ullage) Hg feet high that is quickly pressurized by gaseous 
hydrogen to a constant pressure of 60 pounds per square inch gage. With the gas tem- 
perature held constant at 300° R, the fraction of the total mass of gas that will condense 
and must be replaced in order to hold a constant pressure during the first 5 minutes 
after the start of pressurization is given as follows (based on fig. 7): 


Mass condensed 
Total mass of gas 


At 


0. 17 


60 


H 

^g g 


0. 047 H 


g 


For ullage heights Hg of 0. 5, 1, and 5 feet, the fraction of the total mass that will con- 
dense is 60, 30, and 6 percent, respectively. This estimate indicates that the effect of 
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condensation of the gas should be considered in bulk gas problems such as pressurized 
outflow. 

The effect of interfacial heat transfer from the gas on the bulk gas conditions is now 
considered. For the worst situation described in the preceding paragraph, the inter- 
facial heat transfer from the gas as a fraction of the total enthalpy of the gas, based on 
the analytical data of figure 9 is: 


Interfacial heat f lux from gas 
Total enthalpy of gas 



pi H 
H S g g 




18-5- 

60 


0. 047 x 1020 x 0. 5 


^ 0. 06 or 6 percent 

This conservative estimate indicated that the effect of interfacial heat transfer from the 
gas on the bulk conditions of the gas is small enough to neglect in most bulk gas problems 
and makes the inability of the analysis to predict it of lesser consequence. 

A representative high heat flux from the environment for booster vehicles is about 
50 Btu per square foot per hour. According to figure 8, interfacial heat transfer to the 
liquid is not generally small compared with this figure, so that it should generally be con- 
sidered in bulk liquid heating analyses. 


Effect of Bulk Fluid Mixing on Hydrogen 

The bulk gas temperature used in the analysis T a( j was for a situation where there 
was little gas mixing. It could conceivably be quite different if there were sufficient gas 
mixing, caused by a poorly distributing diffuser, to raise T p above T a( j. Consequently, 
constant values of T f greater than T a( j were used in the analysis to investigate their 
effect on interfacial mass transfer. The result of this study is shown in figure 15. As 
T r increases, there is a tendency toward evaporation, away from the condensation re- 
sult that occurs when the adiabatic temperature is used. With the proper choice of bound- 
ary conditions, no net mass transfer is possible. For the case of no net mass transfer 
(refs. 6 to 8) equations (20) and (29) reduce to 
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(39) 


T r- T sat _ P t /°A 1/2 % 

T sat ” T b p g\“g, / C p g 

This result for zero net mass transfer, which is a function of the fluid properties alone 
(or since the properties are evaluated at T ga j., a function of pressure for a given fluid), 
is plotted in figure 16 to indicate regions of evaporation and condensation. 

The experiment of reference 5 is described here to show that mixing is attainable and 
can have an appreciable effect on interfacial transfer. The results of this experiment 
will be qualitatively compared with figure 16. Stirrers were placed in the liquid and gas 
volumes, near the interface, of a closed Dewar containing two-phase nitrogen. These 
stirrers were used to mix either the liquid or gas volumes in order to determine the ef- 
fect of mixing on the tank pressure. It was observed that the pressure rose when the gas 
was mixed and that the tank pressure dropped when the liquid was mixed. Qualitatively, 
this mixing raises the gas temperature near the interface T r by mixing cool gas near 
the interface with warm gas above it, and lowers the liquid temperature T^ by mixing 
cool liquid of the bulk with warm liquid near the interface. Consider figure 16, which is 
the locus of pressure-temperature conditions for zero mass transfer for a hydrogen sys- 
tem; a similar curve can be generated for nitrogen. Qualitatively, according to figure 16, 
an increase in T r by mixing the gas would tend to cause evaporation and raise the tank 
pressure, while a decrease in T^ by liquid mixing would tend to cause condensation and 
lower the tank pressure. 

It was shown before that condensation could be a major detrimental effect in that it 
increases the mass of pressurant required and heats the liquid. Evaporation at the inter- 
face would be more desirable since it would tend to cool the liquid and decrease pressur- 
ization requirements. According to the previous discussion it can be implied that evapo- 
ration or at least a lower condensation rate can be achieved by mixing the gas by some 
method that does not also mix the liquid. 


CONCLUDING REMARKS 

The condensation that occurs in a tank containing two-phase hydrogen that is rapidly 
pressurized without mixing creates an appreciable transfer of heat to the liquid which 
must be considered in stratification studies. The analysis included herein can be used 
to determine the interfacial heat transfer. The effect, on pressurization requirements, 
of interfacial heat transfer can probably be neglected, while the effect of interfacial mass 
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transfer (condensation) should be considered. Interfacial mass transfer can be deter- 
mined by the analysis contained in this report. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 20, 1965. 
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APPENDIX - SYMBOLS 


A area, sq ft 

Cp specific heat at constant pres- 

sure, Btu/(lb mass)( °R) 

D tank diameter, ft 

H height, ft 

t) specific enthalpy, Btu/lb mass 

h liquid depth (height above aver- 

age tank bottom), ft 

K (Q/A) w /3600p f c p ^D, °R/sec 


Z distance below the interface, ft 

a thermal diffusivity, sqft/sec 

y ratio of specific heats 

rj . dimensionless similarity vari- 

J able, = Z/(4o-.t) 1 / 2 

6 ■ dimensionless temperature, 

_ C p/ T i - T »> 

L 

p density, lb mass/ cu ft 


k 


L 


m/A 


P 

P* 

<Q/A) W 

<V A > 

(q £ A) 


T 

t 


w 


thermal conductivity, Btu/ 

(ft)(°R) or Btu/ 

(ft)(°R)(sec) 

heat of vaporization, Btu/ 
lb mass 

interfacial mass flux, lb mass/ 
(sq ft) (sec) 

tank pressure, psfa 

tank pressure (gage), psig 

heat flux through tank walls from 
environment, Btu/(sq ft)(hr) 

interfacial heat flux from gas, 
Btu/(sq ft)(sec) 

interfacial heat flux to liquid, 
Btu/(sq ft)(sec) 

temperature, °R 

time (after pressurization), sec 

nondimensional velocity, 

= W j (t/ "j )1/2 

flow velocity in Z -direction rela- 
tive to moving interface, ft/sec 


Subscripts: 
ad adiabatic 

b bulk liquid 

g gas 

f final 

int saturation interface in coordinate 

system in h 

i saturation interface in movable 

coordinate system in Z 

j gas or liquid 

£ liquid 

o initial 

r bulk gas 

sat saturation 

sat* helium- hydrogen effective satura- 

tion temperature assumed to 
be that of pure hydrogen 

w wetted wall 

<=° environment outside tank 
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TABLE I. - TESTS PERFORMED 


[Liquid, hydrogen; liquid depth, 55 percent full or 4 ft.] 


Test 

Pressurant 

Type of 

Heat flux 

Nominal pressure- 

Atmos- 


gas 

pressuriza- 

from 

time conditions 

pheric 



tion 

environment, 


pressure 




oa/A) w , 


(nominal), 




Btu/ (sq ft)(hr) 


in. Hg 

1 

Hydrogen 

Artificial 

8 

t < 0: P Q = 0 psig (vented); 

29. 25 

2 

Hydrogen 

Artificial 

8 

t > 0: P* = constant 

29. 24 

3 

Helium 

Artificial 

8 

= 8, 20, 40 psig 

29. 26 

4 


Self 

6.7 

t < 0: P* 0 = 0 (vented); 

28. 78 





t > 0: P* = P’(t) up to 40 psig 




-Diffuser 


rf ^ r ^ 

f y\ '■■- 




SBulk liquids 


-Diffused pressurizing gas 

" — ( q g/ A )i ,nterfacial heat 

transfer from gas 

^"( q l/A)i Interfacial heat 
transfer to liquid 
— (m/A)j Interfacial mass 
transfer 

A-Heat flow 


to side walls 


li . 

l ' 
\ l / 


(Q/ a ) Heat flux from 
v ' w environment 


' / to bottom 

I 1 CD ”8258 

l I L Final equilibrium temperature profile, T sat <P> 

1 L- Nonequilibrium temperature profile, T<Z,t, R> 

'-Initial equilibrium temperature profile, T sat <P 0 > 

Figure 1. - Changes occuring when a tank containing a two-phase fluid is pressurized. 
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Height above tank bottom, 



(a) Tank pressure, 8 pounds per square inch gage. 


Figure 5. - Experimental temperature profile near interface for test 1. 




Height above tank bottom, 




Height above tank bottom, 


50.0 



(c) Tank pressure, 40 pounds per square inch gage. 
Figure 5. - Concluded. 


Height above tank bottom, h, in. 



Temperature, T, °R 


(a) Tank pressure, 8 pounds per square inch gage. 


Figure 6. - Experimental temperature profiles near interface for test 2. 









Interfacial mass flux (condensation), (m/A)j, lb/(sq ft)(hr) 







Height above tank bottom, 



(a) Tank pressure, 8 pounds per square inch gage. 

Figure 10. - Experimental temperature profiles near interface for test 3 with helium pressurant. 
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rture, T, “R 

pounds per square inch gage. 
). - Concluded. 



Height above tank bottom, h, in. 


52 



Figure 11. - Experimental temperature profiles near interface for test 4 with self-pressurization. 
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Interfacial mass flux (condensation) (m/A)j, lb/(sq ft)(hr) 
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0 . 5 . 
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Ratio of gas to liquid temperature gradients at 
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Figure 16. - Direction of mass flux as function of bulk fluid tempera- 
tures and tank pressure for hydrogen. 
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